Stripe formation in differentially forced binary systems 
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We consider pattern formation in periodically forced binary systems. In particular we focus on 
systems in which the two species are differentially forced, one being accelerated with respect to 
the other. Using a continuum model consisting of two isothermal ideal gases which interact via a 
frictional force we demonstrate analytically that stripes form spontaneously above a critical forcing 
amplitude. The wavelength of the stripes is found to be close to the wavelength of sound in the limit 
of small viscosity. The results are confirmed numerically. We suggest that the same mechanism may 
contribute to the formation of stripes in experiments on horizontally oscillated granular mixtures. 



Binary systems subject to an oscillatory driving force 
are often found to phase segregate and form patterned 
structures 0, Q . Often this segregation can have impor- 
tant practical applications. For example it may provide a 
useful way of separating two mixed species, or conversely, 
it may be undesirable in an industrial process 0. The 
phenomenon is well documented experimentally but a full 
theoretical understanding is still missing with different 
mechanisms for the pattern formation being proposed in 
the literature. 

Sanchez et al. 4] have carried out experiments and 
parallel simulations on granular mixtures immersed in 
water. For sufficiently large amplitudes of vibration a 
mixture of equal sized glass and bronze particles were 
found to separate into a striped pattern. Sanchez et al. 
proposed that the differential fluid drag on the two com- 
ponents of the mixture is the mechanism responsible for 
the segregation. 

Mullin et al. 0, 0, 0| have performed a series of exper- 
iments in which they horizontally oscillate a quasi-2 di- 
mensional layer of bronze spheres and poppy seeds. The 
species are again observed to segregate into a striped pat- 
tern. However they propose a different explanation for 
the stripes: the effective excluded volume interactions 
which occur because of the different size of the shaken 
particles 0. 

Therefore our aim in this Letter is to investigate possi- 
ble mechanisms for binary phase segregation by providing 
an approximate analytic solution to a one-dimensional 
isothermal continuum model describing the physics of a 
binary mixture of two ideal fluids. We find that if the 
components are differentially forced stripes are formed 
in the concentration above a critical forcing amplitude. 
The wavelength of the stripes is controlled by the veloc- 
ity of sound and standing wave oscillations in the total 
density pay a crucial role in the stripe formation. 

We then relate our results more closely to experiments 
on granular mixtures through a simple particle model. 
We find that differential forcing of the two types of par- 
ticle is enough to cause stripe formation even if the par- 
ticles have the same mass and volume. 

We consider two ideal fluids, labeled A and B, which 



are coupled by a frictional force proportional to the differ- 
ence in their velocities. The physical origin of this force 
is from collisions between the A and B particles which 
tend to equalize the velocities. The A fluid is acceler- 
ated by a periodic force acos(wi). The Navier-Stokes 
and continuity equations for this system may be written 



p D t u 



u A )+d x (vp A d x u A ) 



6d xP A + £(u B 
+ap A cos(wi), (1) 
p B Dfu B = -6d xP B + £(u A - u B ) + d x (vp B d x u B ),(2) 
Dfp A = -p A d x u A , (3) 
D B p B = -p B d x u B , (4) 



where — ksT/m is the reduced temperature, v is the 
kinematic viscosity, and p A , p B and u A , u B are the densi- 
ties and velocities of the A and B fluids. D A = dt + u A d x 
and D B = dt + u B d x are the material derivatives in the A 
and B frames, and the friction coefficient can be written 
£ = £° p A p B I p, where £° is constant and p = p A + p B is 
the total density of the mixture. 

Fig. H shows the results of numerically solving 
Eqns. (^QJ using a lattice Boltzmann algorithm. Ini- 
tially the densities of the two components have the same 
constant value po plus a small random perturbation. Pro- 
vided the amplitude of the forcing a is sufficiently high, 
strongly segregated A and B stripes form as shown in Fig. 

These correspond to an approximately sinusoidal 
modulation of the concentration <f> = p A /p. We define 
the average velocity of the fluid u = (p A u A + p B u B )/p 
and the rest frame of the fluid to be that moving at ve- 
locity (u), where (...) denotes a spatial average. In this 
frame the total density performs standing wave oscilla- 
tions at the forcing frequency but shifted in phase. 

We assume that the coupling between the two com- 
ponents is sufficiently strong that it is the fastest relax- 
ation process in the system. Therefore u A and u B are 
always very close in value; hence their time derivatives 
are essentially the same. This allows us to approximate 
D A u A = D B u B = D t u, where D t = d t + ud x . Adding 
equations and (J2J to remove the coupling term then 
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FIG. 1: Evolution in time of (a) the concentration and (b) 
the total density of a differentially forced binary fluid as a 
function of the distance x along a one dimensional system, 
observed in the rest frame of the fluid. The concentration 
forms stripes and the total density performs standing wave 
oscillations, which increase in amplitude, at the frequency of 
the applied force. The parameters used were £° = 6, v — 0.5, 
= 1, a = 0.358, and ui — 0.2. The wavelength of the 
resulting stripes A = 33 is close to the wavelength of sound 



gives 



31. 



D t u — — d x p + a<f> cos( cut) + vd 2 u. 
P 

Similarly adding Eqns. JSJ and (0J 

D t p = -pd x u. 



(5) 



(6) 



Motivated by the numerical results we consider the trial 
solutions 



-+A0sin(foc), 



(7) 



p = 2po + cos(kx) [Api cos(wi) + Ap2 sm(u>t)] , (8) 
a sin(wf) 
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-sm(kx) [Aui cos(wt) + Au 2 sin(wi)] ,(9) 



where k = ^ is the wavenumber of the stripes. 

Substitution of the trial solutions (QQ into Eqns. (JSJ 
and JJjJ, and comparing cosine and sine terms gives 
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The crucial information needed is whether a wave with 
wavenumber k grows or decays, and the rate of this pro- 
cess. To understand this we must look at the behaviour of 
<p as a function of time. The material derivative must now 
be treated exactly because the concentration changes 
much more slowly than densities and velocities. Using 
the continuity equations 10 and and p A as p B p 
gives 



D, 



1 



d x (u B ~u A ) + 



-d x p. (14) 



4 0,00 
Rearranging Eqn. J2J and substituting the trial solutions 
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Apt ) cos(urf). (16) 



Comparing the first order terms in ub — ua in Eqn. 115|) 
and in u in Eqn. © we see that the condition that A and 
B are strongly coupled is equivalent to stating £° >> lu. 
Substituting © and (JTSJl into JT3I) leads to 

if Ck 
D t <j> = — -9A4>k 2 sin(fcx) + — cos(fca;) 
s V 2 

+ J— - coa(uit) sin(fcx) [Api cos(oji) + Ap2 sin(ujt)] 
opo 

The time scale over which <f> changes is much longer than 
that taken for a single oscillation of the force. This allows 
us to average the right hand side over time 
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Substituting expression (jlt)ft for Api and integrating 
gives the exponential solution 



A<t> 



(18) 
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FIG. 2: (a) Growth rate of stripes T as a function of wavenum- 
ber k (x simulations, theory), (b) Critical forcing ampli- 
tude a c , above which the system forms stripes, as a function 
of angular frequency ui for two different viscosities ( x simula- 
tions with v = 2, o simulations with v = 0.5, ■ ■ • theory with 

v = 2, theory with v = 0.5). (c) Stripe wavelength A as 

a function of speed of sound c s for an angular frequency of 
uj — 0.4 and different values of viscosity v ( — wavelength of 
sound) . 
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A typical curve of T against k is shown in Fig. E^a) . 
The peak represents the fastest growing wavenumber. 
The numerical results, presented for comparison, were 
obtained by solving Eqns. Q-lQJ for the parameters 



At = 0.001, uj = 0.2, £° = 32, v = 2, 6 = f, p Q = 4, 
Acj) — 0.001, a — 0.8 and k = 2n/L such that one stripe 
filled the entire system of length L. The initial condition 
was taken to be the trial solution. 

Using (QU we may make two predictions. First, stripe 
formation only occurs when the amplitude of the force is 
above a certain critical value a c such that T is positive 
for some value of k. To test this numerically we used 
a system with L = 1500, At = 0.01, £° = 6 and = 
1. Results are shown in Fig. Etb). In each simulation 
we initially set the density to be po = 4 plus a small 
random perturbation of amplitude 0.02. Stripe formation 
was defined to occur when ((<fi ~ 0o) 2 ) > 0.01 and the 
system was observed for up to 20000 time steps. Exact 
agreement between the analytic and numerical results is 
not expected because the simulation size is finite and the 
value of k is limited to integer multiples of • Moreover 
the assumption that £° >> ui breaks down for large u>. 

Second, the wavelength of the stripes tends toward the 
wavelength of sound A s = 2ttc s /uj, where c s = yfd for an 
isothermal gas, in the low viscosity limit. This is clearly 
illustrated in Fig.|2Jc), which show how stripe wavelength 
changes as a function of the speed of sound for systems at 
a c . The relation to c s is expected because the evolution 
of the total density ijBJ) is a standing wave (i.e. a sound 
wave) which slowly increases in amplitude as the striped 
concentration profile forms. 

Thus we have shown that, within a continuum theory, 
differential forcing of the two components of a binary 
fluid can lead to stripe formation in the concentration 
profile. 

We now take the first steps in investigating whether 
the same mechanism can hold in granular systems by so- 
lution of a two-dimensional, particle-based model which 
includes an excluded volume constraint. We consider two 
types of particles A and B which have equal diameter d, 
and equal mass m. The particles undergo free streaming 
and a periodic force is applied to the A particles. 
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(20) 
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Each particle is then identified as belonging to one of 
a grid of square cells with dimension 2ri. Each particle 
is checked to see if it overlaps with particles in its own 
cell or the three other closest cells. If two overlapping 
particles are travelling towards each other then they col- 
lide elastically, thus conserving energy and momentum. 
The smaller the time-step At the more rigorously the ex- 
cluded volume condition is enforced. The system is made 
isothermal by dividing it into slices of width 2ri perpen- 
dicular to the forcing direction and rescaling velocities 
in the center of mass frame of each slice to accord with 
equipartition of energy. 

The system forms stripes of A and B particles per- 
pendicular to the forcing direction as in the continuum 



4 



' ■ ■ 



mm 



■f^>'< •••§;.>v"--- . . - ■ ■ "' 

... -i . J-.-:-.v , 



(a) t=0.25s 



(b) t=2.5s 



(c) t=25s 



FIG. 3: Time evolution in seconds of a horizontally shaken 
granular system. The different friction coefficients of the par- 
ticle species cause the stripe formation. The parameters of 
the simulations are given in the text. 



model. Moving to two dimensions and including fluctu- 
ations does not qualitatively alter the stripe formation 
mechanism. However an interesting difference which we 
do not yet have an explanation for is that there is no 
threshold forcing amplitude for the stripe formation. In- 
stead the time taken to form stripes diverges as the forc- 
ing tends to zero. The wavelengths of the stripes corre- 
sponds to the wavelength of sound to within ~ 10%. 

Finally we discuss how differential forcing of the two 
components of a granular mixture may be responsible for 
the experimental findings of Mullin et al. |5j . Mullin et 
al. consider a flat tray which is oscillated in the horizon- 
tal plane at a frequency of \2Hz. A mixture of poppy 
seeds and bronze spheres are placed on the tray. The 
particles are initially mixed and form an almost com- 
plete single particle layer. Forcing of individual particles 
arises principally from the frictional interaction with the 
base of the container. Mullin et al. observe that in time 
the two types of particles segregate into a striped pattern 
(very similar to that in Fig. |3 which shows the results of 
the computer simulations described below). 

We modelled this system using the excluded volume 
elastic scattering algorithm with the update in the par- 
ticle velocities given by 



A.B 



(t+At) 



A,B 



(t) + (acos(ujt)- f j, A - B gv l ) At, (22) 



where a is the amplitude of the acceleration applied to 
both species. The final term represents the frictional in- 
teraction with the container base: g — 9.8ms~ 2 is the 
acceleration due to gravity, Vj is the unit vector in the 
direction of the particle velocity, and /i A ' B are the co- 
efficients of friction for poppy seeds and bronze spheres 
respectively. No thermostat is used. Energy from the 
forcing is dissipated by friction. Because the friction co- 
efficients of the two types of particles are different they 
are differentially forced, although this difference at large 



a is like a square wave rather than the sinusoidal forcing 
used in Eqn. J3TJ . 

We make the approximation that the static and dy- 
namic friction coefficients are equal and ignore questions 
as to how much the spheres roll or slide. Experimental 
data for /u, A and fi B is not available and we use the phys- 
ically realistic values [i A — 0.6 and [i B = 0.2. We find 
that the qualitative behaviour is not dependent on the 
exact values. 

Fig. |3 is the result of a simulation showing the evolu- 
tion in time of a shaken binary system which is initially 
mixed. The dimensions of the tray are 90mm by 180mm 
and there are 1500 particles each with radius r = 1.5mm. 
The oscillation frequency was 4.8Hz and the forcing am- 
plitude a — g. Stripes form on a time scale very similar 
to that found in the experiments, which was ~ 40s 
Therefore it seems reasonable that the differential ac- 
celeration of the two components can be identified as at 
least partially contributing to the stripe formation. How- 
ever we caution that these simulations are not isothermal, 
hence direct comparison with the isothermal model pre- 
sented earlier is difficult. Furthermore the model assumes 
elastic collisions which may not be a good approximation 
for the real experiment. 

To summarise, we have demonstrated analytically that 
if two isothermal ideal gases interacting via a frictional 
force are subject to differential forcing stripes form in 
the concentration profile above a critical forcing ampli- 
tude. The stripe formation originates from energy be- 
ing pumped into standing wave oscillations, and conse- 
quently the wavelength of the stripes is approximately 
that of the wavelength of sound. These calculations were 
carried out for a continuum model. We then used a 
particle-based numerical model to provide evidence that 
the same mechanism, differential forcing due to differen- 
tial friction with the base of the container, is likely to be 
responsible for stripe formation in horizontally oscillated 
granular mixtures. Further work will aim to use more 
realistic models of granular systems to investigate the ef- 
fect of inelastic collisions and non-uniform temperature 
distributions on the pattern formation. 

We thank M.R. Swift, T. Mullin and C. Dcnniston for 
helpful discussions. 
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